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Summary  of  results 

The  target  of  the  work  is  to  use  Transformation  Optics  (TO)  in  order  to  guide  as  well  as  focus  plasmons.  In  the 
initial  phase  of  the  work  we  focused  on  using  TO  to  GRaded  INdex  (GRIN)  Lenses  of  Luneburg  type  and  explored 
focusing  and  guiding.  In  particular  we  showed  both  analytically  and  numerically  with  home-made  codes  that  for  plane 
wave  initial  conditions  (a)  A  sequence  of  Luneburg  lenses  can  form  a  waveguide  with  low  losses  that  can  guide  EM 
waves  in  linear  as  well  as  curved  geometries,  (b)  Random  distributions  of  Luneburg  lenses  can  lead  to  the  formation 
of  caustics  and  optical  rogue  waves,  (c)  We  used  commercial  code  GOMSOL  and  generated  plasmon  waves  at  a 
metal-dielectric  interface.  Subsequently  we  analysed  propagation  of  Gaussian  beams  and  found  that  Luneburg  lens 
focusing  is  not  as  good  as  in  the  case  of  plane  waves.  In  this  case,  the  addition  of  nonlinearity  of  the  Kerr  type 
improves  focusing  and  thus  propagation  in  the  Luneburg  waveguides.  In  the  case  of  plasmon  propagation  we  found 
that  Luneburg  lenses  may  focus  plasmons  in  a  resonant  fashion. 

The  results  of  the  Luneburg  work  formed  a  manuscript  that  has  been  invited  to  a  special  volume  to  be  published 
in  the  Journal  of  Optics.  This  article  is  currently  under  review.  We  anticipate  the  preparation  of  one  more  article 
with  the  nonlinear  effects. 
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I.  INTRODUCTION 

Gradient  Index  (GRIN)  metamaterials  are  formed  through  spatial  variation  of  the  index  of  refraction  and  lead  to 
enhanced  light  manipulation  in  a  variety  of  circumstances.  These  metamaterials  provide  natural  means  for  construct¬ 
ing  various  types  of  waveguides  and  other  optical  configurations  that  guide  and  focus  light  in  specific  desired  paths. 
Different  configurations  have  been  tested  experimentally  while  the  typical  theoretical  approach  uses  Transformation 
Optics  (TO)  methods  to  cast  the  original  inhomogeneous  index  problem  to  an  equivalent  one  in  a  deformed  space 
[2]  [3]  [4]  [5]  [6] .  While  this  approach  is  mathematically  elegant,  it  occasionally  hides  the  intuition  obtainned  through 
more  direct  means.  Furthermore,  a  general,  continuous  GRIN  waveguide  may  be  hard  to  analyse  in  more  elemental 
units  and  relate  its  global  features  to  these  units.  In  the  present  work,  we  adopt  precisely  this  latter  avenue,  viz. 
attempt  to  construct  waveguide  structures  that  are  seen  as  lattices,  or  networks,  of  units  with  specific  features.  This 
is  a  ”  metamaterials  approach” ,  where  specific  properties  of  the  ”  atomistic”  units  are  inherited  as  well  as  expanded  in 
the  network. 

The  ’’atomic”  unit  of  the  network  is  a  Luneburg  lens  (LL);  the  later  is  a  spherical  construction  where  the  index  or 
refraction  varies  from  one,  in  its  outer  boundary,  to  \f2  in  the  centre  with  a  specific  functional  dependence  on  the  lens 
radius  [7].  Its  basic  property,  in  the  geometrical  optics  limit,  is  to  focus  parallel  rays  on  the  spherical  surface  on  the 
opposite  side  of  the  lens  [7]  [8]  [9]  [10].  This  feature  makes  LL’s  quite  interesting  for  applications  since  the  focal  surface 
is  predefined  for  parallel  rays  of  any  initial  angle.  While  the  rays  traverse  the  lens,  they  suffer  variable  deflections 
depending  on  their  distance  from  the  optical  axis  ray,  leading  to  a  point  image  on  the  lens  surface.  This  property 
of  the  spherical  LL  is  also  shared  by  its  cylindrical  equivalent  formed  by  long  dielectric  cylinders  while  the  light 
wavevector  impinges  perpendicularly  to  the  cylindrical  axis.  This  geometry  turns  the  problem  into  a  two  dimensional 
one,  constructed  for  any  plane  that  cuts  the  LL  cylinder  perpendicular  to  its  axis.  The  work  in  this  article  focuses  on 
exactly  this  type  of  cylindrical  LLs  and,  as  a  result,  our  approach  is  strictly  two  dimensional  [7]  [8]  [10]. 


II.  SINGLE  LUNEBURG  LENS 
A.  Ray  tracing  solution 

The  radial  dependence  of  the  refractive  index  n  of  a  planar  LL  is  given  by  [7] 


n(r)  =  y2-(^)2  (1) 

where  R  is  the  radius  of  lens  while  r  is  the  radial  distance  from  the  center  in  the  interior  of  the  lens.  In  order  to 
find  an  analytical  solution  for  the  ray  path  inside  a  LL  we  use  Fermat’s  principle  for  optical  path  minimization  of  the 
action  S  and  find  [1] 


r(0) 


_ qi? _ 

^1-  sin(2((/.  + /?)) 


(2) 


where  C  and  f3  are  constants.  Using  appropriate  boundary  solutions  we  may  obtain  analytically  the  ray  propagation 
in  a  single  Luneburg  lens.  In  Fig.  1  we  plot  several  rays  through  the  analytical  solution;  it  is  clear  that  the  solution 
describes  fully  the  physics  of  the  rays  inside  a  single  LL. 


Distribution  A:  Approved  for  pubiic  reiease;  distribution  is  uniimited. 


4 


single  Luneburg  lens,  Initial  propagation  angle  0=0 


single  Luneburg  lens,  Initial  propagation  angle  0=7t/4 


FIG.  1:  A  Luneburg  lens  is  shown  with  the  red  dashed  line.  The  blue  lines  depict  parallel  rays  for  initial  propagation  angles 
(i)  ^  =  0  (upper  figure)  and  (ii)  ^  =  7r/4  (lower  figure).  All  rays  are  focused  onto  a  single  point  on  the  opposite  side  of  the 
Luneburg  lens  leading  to  thus  we  have  perfect  imaging. 


III.  LUNEBURG  LENS  WAVEGUIDES 

The  analytical  solution  for  the  ray  propagation  through  a  single  LL  may  be  used  in  order  to  study  analytically 
the  ray  transfer  through  various  configurations  of  LLs  that  form  waveguides.  Using  the  initial  entry  point  on  the  LL 
(xo,^o)  as  well  as  the  initial  ray  angle  0  the  exit  point  {x^y)  and  the  associated  exit  angle  O'.  We  may  thus  form  a 
mapping  from  {xo^yo^O)  to  {x^y^O');  further  propagation  in  the  surrounding  medium  is  rectilinear  while  the  entry  to 
the  next  LL  is  governed  by  a  new  initial  entry  point  with  angle  equal  to  the  previous  exit  angle.  The  resulting  ray 
may  be  traced  quite  easily. 

The  first  arrangement  we  use  a  series  of  touching  LLs  on  a  straight  line  as  in  Fig.  2;  we  note  that  due  to  time 
reversal  symmetry  the  focusing  point  of  a  LL  produces  a  ray  bunch  that  exits  mutually  parallel  form  the  next  one. 
As  a  result  we  have  the  formation  of  a  waveguide,  that,  depending  on  wheather  the  number  of  lenses  is  even  or  odd, 
defocusses  or  focusses  respectively. 

The  analytical  solution  may  be  used  for  curvilinear  arrangements  of  LLs  as  well.  In  Fig.  3  we  show  a  network  of 
21  Luneburg  lenses  that  are  placed  in  such  a  way  so  that  to  bend  light  almost  at  right  angle.  While  there  are  some 
ray  losses  in  the  bend,  the  vast  majority  of  the  rays  pass  through  leading  to  an  efficient  curved  waveguide.  Choosing 
an  odd  total  number  of  lenses  leads  to  focussing  in  the  last  lens,  while  an  even  number  of  LLs  would  give  rise  to 
defocussing,  as  in  the  linear  case. 
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Linear  Luneburg  waveguide,  formed  by  7  lenses 


Linear  Luneburg  waveguide,  formed  by  6  lenses 


FIG.  2:  The  red  dotted  lines  denote  the  arrangement  of  lenses.  The  blue  lines  show  the  ray  tracing.  Light  is  guided  by  Luneburg 
lenses  accross  the  linear  network,  (i)  Arrangement  with  7  lenses  (upper  figure).  We  can  see  that  all  rays  are  focus  on  the  last 
lens,  (ii)  Arrangement  with  6  lenses  (lower  figure).  All  rays  exit  in  the  same  mode  as  originally  entered,  i.e.  parallel  to  the 
waveguid  axis. 


Bending  a  beam,  arrangement  of  21  lenses 


FIG.  3:  The  red  dotted  lines  denote  the  network  of  Luneburg  lenses,  while  the  blue  lines  show  the  rays.  Except  for  some 
rays  that  escape  from  the  waveguide,  the  majority  of  rays  is  guided  through  the  arrangement  and  exit  at  right  angles  wrt  the 
oribinal  direction.  The  rays  are  focused  in  the  last  lens  showing  the  also  in  curved  geometries  waveguides  with  odd  number  of 
lenses  focus  the  rays  while  waveguides  with  even  numbers  of  LLs  do  not. 


IV.  PARAMETRIC  RAY  SOLUTION  FOR  SINGLE  LUNEBURG  LENS 

The  approach  we  used  [1]  is  based  on  an  equation  that  is  not  a  real  dynamical  equation  for  a  two  dimensional 
(2D)  problem,  since  we  have  chosen  as  generalized  time  one  of  the  space  coordinates;  i.e.  the  radial  variable  r.  The 
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advantage  of  this  choice  is  that  it  leads  to  a  one  dimensional  Lagrangian,  that  results  in  equation  of  motion  much 
easier  than  if  we  had  a  true  2D  Lagrangian  [1].  Unfortunately,  the  solution  of  this  equation  being  a  function  of  ordinate 
as  a  function  of  abscissa  fails  to  work  when,  for  a  particular  ray  or  LL  arrangement,  the  ray  must  bent  backwards. 
The  solution  to  this  problem  is  to  find  a  true  2D  equation  of  motion  by  using  generalized  time  in  addition  to  the  two 
space  coordinates;  the  Lagrangian  in  this  case  contains  two  spatial  coordinates.  Unfortunately,  the  resulting  equations 
are  quite  complex  and  we  thus  use  numerics  for  their  direct  solution.  It  is  preferable  to  use  Hamilton’s  equations  in 
Cartesian  coordinates;  the  latter  are  first  order  and  thus  easier  to  handle  numerically,  the  specific  expressions  can  be 
found  in  ref.  [1];  we  show  here  the  associated  figure  for  a  single  Luneburg  lens. 


single  Luneburg  lens,  initial  propagation  angle  0=7t/3 


FIG.  4:  Use  of  the  Hamiltonian  ray-tracing  method  and  through  numerical  solution  of  the  associated  diferential  equations  we 
obtain  the  propagation  of  light  through  a  single  Luneburg  lens  with  initial  propagation  angle  0  =  7r/5. 


V.  WAVE  PROPAGATION  IN  LUNEBURG  LENS  NETWORKS 

We  may  compare  the  analytical  ray  propagation  results  to  simulations  that  solve  directly  Maxwell  equations.  We 
use  the  finite  difference  in  time  domain  (FDTD)  method  based  on  Yee’s  algorithm  developed  by  Allen  Tafiove  [13]. 
Specifically,  FDTD  handles  light  as  an  electromagnetic  field  since  it  solves  numerically  the  fully  time  dependent 
Maxwell  field  equations: 


dE 

dt 

dH 

dt 


-V  X  # 

xE 


(3) 

(4) 


where  E  and  H  denote  the  electric  and  magnetic  fields  respectively,  while  e  and  /i  are  the  dielectric  permittivity  and 
the  magnetic  permeability  of  the  medium.  Specifically,  we  consider  that  the  material  outside  of  the  LL  lens  is  air, 
thus  setting  5  =  5o  and  /i  =  /io,  where  we  have  normalized  5o  =  /io  =  1-  Additionally,  inside  the  lens  the  permeability 
is  maintained  the  same,  since  we  consider  that  the  lenses  are  fabricated  by  optical  materials.  On  the  other  hand 
dielectric  permittivity  is  a  function  of  the  position  and  we  can  obtain  it  from  the  square  of  the  Luneburg  refractive 
index,  viz.  Eqn.  (1),  thus 


“  ^  R2  -  2 


R2 


We  consider  TM  modes  propagating  in  the  x'x  axis,  i.e. 


Ez  —  Ey  —  Hx  —  0 


(5) 

(6) 
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As  a  result,  according  the  condition  of  the  Eq.  (6),  the  Eqns.  (3)  and  (4)  take  the  following  form 

dHy  1  dEx 

dt  fi  dz 

dHz  1  dEx 

dt  /i  dy 

dEx  ^  1.5^  _  dHy 
dt  s  dy  dz 

Eurthermore  we  use  a  source  creating  plane  waves  [13].  The  source  is  located  in  the  front  of  the  lattice,  i.e. 

E(l,  j)  =  Eosin{Lut) 


(7) 

(8) 
(9) 


(10) 


where  the  index  (l,j)  denotes  that  the  source  is  located  on  the  first  section  of  the  lattice  in  the  y-direction,  uj  is  the 
angular  frequency  of  the  EM  wave  and  it  is  defined  as  a;  =  2'Kcj\  while  Eq  is  the  amplitude  of  the  wave. 

In  the  lattice  edge  we  apply  absorbing  boundary  conditions  (ABCs)  [13],  while  we  perform  all  computations  in  the 
microwave  regime.  Specifically,  we  take  the  wavelength  of  the  EM  wave  as  A  =  1cm.  Einally,  the  radius  of  each  lens 
in  the  following  simulations,  is  six  times  greater  that  the  wavelength,  thus  R  =  6X  =  6cm. 

In  Eigs.  (5-8)  we  show  the  results  for  the  full  propagation  of  EM  waves  through  LLs  and  LLWs.  Eor  each  figure, 
in  the  upper  diagram  we  plot  the  steady  state  of  the  intensity  of  the  electric  field,  while  in  the  lower  diagram  we  plot 
the  electric  field  for  a  specific  time,  i.e.  we  take  a  snapshot  of  the  electric  field.  In  order  to  reach  steady  state  of 
the  intensity,  we  integrate  the  square  of  the  electric  field  for  one  period.  Additionally,  in  order  to  be  sure  that  the 
integration  is  made  after  the  EM  wave  has  passed  through  the  Luneburg  network,  we  choose  to  average  after  long 
time,  i.e. 


/  =  lim 


nt-\-27T 

Jt 


E{r,  t)^dt 


(11) 


In  Eig.  (5)  we  show  the  results  for  a  single  LL  where  we  observe  focussing  on  the  lens  surface  in  agreement  with 
both  the  theoretical  prediction  of  Eig.  (I)  and  with  Hamiltonian  ray  tracing  method  of  Eig.  (4). 
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Single  Luneburg  lens,  R=6?l,  screenshot  of  Electric  Field 


FIG.  5:  Simulation  with  FDTD  method  for  a  single  Luneburg  lens.  The  EM  wave  is  focused  on  the  surface  of  the  lens,  (i) 
(Upper  figure)  Plot  of  the  steady  state  of  the  intensity  /  of  the  electric  field,  (ii)  (Lower  figure)  Plot  of  the  electric  field  for  a 
specific  time,  i.e.  a  snapshot  of  the  electric  field  E. 


Subsequently  we  ivestiagate  the  linear  arrangement  of  LL  waveguids  with  six  and  seven  lenses  respectively  (Figs. 
(2)  and  (5)).  In  case  with  even  number  of  lenses,  as  in  Fig.  (6),  the  ray  beam  is  guided  and  it  defocuses  in  the  last 
lens.  On  the  other  hand,  in  the  case  with  odd  number  of  lenses,  as  in  Fig.  (7),  the  beam  is  focused  on  the  last  lens. 
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FIG.  6:  Simulation  with  FDTD  method.  Linear  waveguide  arrangement  formed  with  six  Luneburg  lenses.  The  EM  wave  is 
guided  through  the  linear  network  of  lenses  and  the  wave  defocuses  in  the  last  lens.  This  result  is  in  agreement  with  figure  2. 
(i)  (Upper  figure)  Plot  the  steady  state  of  the  intensity  /  of  the  electric  field,  (ii)  (Lower  figure)  Plot  the  electric  field  for  a 
specific  time,  i.e.  a  snapshot  of  the  electric  field  E. 
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Linear  Luneburg  waveguide  formed  by  7  lenses  with  R=6^ 
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Linear  Luneburg  waveguide  formed  by  7  lenses  with  R=6^ 


FIG.  7:  Simulation  with  FDTD  method.  Linear  Luneburg  waveguide  arrangement  with  sevenlenses.  Beam  guiding  and  focusing 
on  the  surface  of  the  last  lens  is  observed.  The  result  is  in  agreement  figure  2.  (i)  (Upper  figure)  Plot  the  steady  state  of  the 
intensity  /  of  the  electric  field,  (ii)  (Lower  figure)  Plot  the  electric  field  for  a  specific  time,  i.e.  a  snapshot  of  the  electric  field 
E. 


In  the  Fig.  (8)  we  investigate  right  angle  light  bending  while  in  Fig. 
We  observe  propagation  and  guiding  with  small  losses. 


(9)  we  have  a  circular  arrangement  of  LL’s. 
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FIG.  8:  A  rectangular  waveguide  is  formed  with  eleven  Luneburg  lenses.  The  EM  wave  is  guided  and  bends  at  right  angles 
through  a  network  of  lenses,  (i)  (Upper  figure)  Plot  the  steady  state  of  the  intensity  I  of  the  electric  field,  (ii)  (Lower  figure) 
Plot  the  electric  field  for  a  specific  time,  i.e.  a  snapshot  of  the  electric  field  E. 
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FIG.  9:  A  circular  waveguide  formed  by  seventeen  Luneburg  lenses.  The  EM  wave  is  guided  through  a  circular  orbit,  (i) 
(Upper  figure)  Plot  the  steady  state  of  the  intensity  I  of  the  electric  field,  (ii)  (Lower  figure)  Plot  the  electric  field  for  a  specific 
time,  i.e.  a  snapshot  of  the  electric  field  E. 


VI.  LUNEBURG  NONLINEAR  EFFECTS 

In  order  to  analyze  specific,  realist c  geometries  we  purchased  and  implemented  in  teh  context  of  the  present  grant  the 
commercial  package  COMSOL.  During  the  initial  learing  phase  we  studied  simple  dielectric  structures.  Subsequently 
we  entered  the  study  of  specific  LL  configurations  as  well  as  investigated  the  role  of  nonlinearity  in  GRIN  media 
made  of  Luneburg  lenses.  While  plane  waves  focus  perfectly  in  LLs,  more  realistic  wavepackets  have  distortions. 
Specifically,  we  tested  initially  gaussian  as  well  as  square  wave  sources  and  found  that  LL  focusing  takes  place  slightly 
outside  the  lens.  This  distortion  can  be  corrected  through  nonlinearity.  The  results  are  shown  in  Fig.  10. 
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FIG.  10:  (a)  Square  wave-single  frequency  pulse  focusing  in  a  single  Luneburg  lens,  (b)  Same  as  in  (a)  but  with  nolinearity 
X  =  1.0  X  10~^^.  We  observe  improvement  in  focusing  induced  through  nonlinearity 

We  may  also  form  similarly  an  LL  waveguide  and  observe  the  imporvement  in  focusing  and  propagation  as  in  Fig. 
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Surface:  emw.normE''2  (m^*kg^/(s**A^)) 
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Surface:  emw.normE''2  (m^*kg^/(s**A^)) 
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FIG.  11: 


(a)  Square  wave-single  frequency  pulse  propagation  in  a  linear  four  Luneburg  lens  waveguide,  (b)  Same  as  in  (a)  but 


with  nolinearity  y  =4^0  x  10 
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VII.  LUNEBURG  PLASMON  FOCUSING 

Finally,  we  used  the  commercial  package  COMSOL  in  order  to  study  the  effect  of  LL  focusing  of  plasmons.  The 
results  show  that  evanescent  field  focusing  results  in  plasmon  wave  concentration.  The  geometry  we  use  involves  two 
half  planes,  the  lower  one  being  a  metal  (glod)  while  the  uper  one  a  dielectric  (air).  In  the  upper  plane  we  place  (half) 
Luneburg  lenses.  What  we  are  interested  in  is  the  focusing  and  manipulation  of  plasnoms  generated  on  the  interface 
through  LLs.  In  Fig  (12)  we  show  wave  focusing  in  the  dielectric  that  results  in  focusing  of  the  surface  plasmon  that 
is  generated  at  the  interface.  IN  Fig  (13)  we  have  an  array  of  seven  lenses  in  the  dielectric  plane.  The  impinging 
microwave  radiation  generates  surface  plasmons  that  propagate  across  the  interface  together  with  the  wave  in  the 
dielectric  part  that  is  being  focused  by  the  LL’s.  IN  Fig.  (14)  we  show  the  wave  amplitude  just  above  and  just  below 
the  interface  accross  the  propagation  direction.  We  observe  the  focusing  effect  of  the  Luneburg  lenses  (a)  while,  at 
the  same  time,  we  see  a  resonant  enhancement  of  the  plasmon  wave  in  the  metal.  This  is  a  very  interesting  effect  that 
we  believe  stems  from  the  matching  of  the  propagating  surface  plasmon  to  the  focusing  of  the  wave  in  the  dielectric 
by  the  Luneburg  lenses.  We  are  currently  studying  this  peculiar  effect  analytically  and  we  believe  that  shortly  we 
will  be  able  to  address  it  quantitatively.  It  is  obvious  that  this  effect  is  potentially  very  interesting  for  plasmon  decay 
supression  during  propagation. 
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FIG.  12:  (a)  Square  wave-single  frequency  pulse  propagation  in  a  single  Luneburg  lens  in  touch  with  a  metal  (gold).  We  observe 
plasmon  focusing,  (b)  Zoom  of  (a)  in  the  focusing  region 


Distribution  A:  Approved  for  pubiic  reiease;  distribution  is  uniimited. 


17 


Surface:  Electric  field  norm  CV/m) 
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FIG.  13:  Generation,  propagation  and  focusing  of  plasmons  in  a  Luneburg  lens  waveguide. 
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Line  Graph;  Electric  field  norm  (V/m) 


Line  Graph;  emw.normE''!  (V/m) 


FIG.  14:  (a)  Electrif  field  distribution  along  the  propagation  direction  of  Fig.  (13)  in  the  dielectric  just  above  the  interfase.  We 
observe  the  Luneburg  induced  focusing,  (b)  Electric  field  plasmon  distribution  just  below  the  inteface.  We  observe  plasmon 
resonant  enhancement. 
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